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Observational evidence for the presence of PAHs in distant 
Luminous Infrared Galaxies using ISO and Spitzer 
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Abstract. We present ISOCAM 15yum and MIPS 24 /jm photometry of a sample of 16 distant Luminous Infrared Galaxies 
(LIRGs) characterized by a median luminosity Lir ~ 2x10" L Q and redshift z = 0.7 (distributed from z = 0.1 to 1.2). 
While some sources display 24/15 yum flux ratios also consistent with a featureless continuum dominating their mid-infrared 
(MIR) spectral energy distributions (SEDs), the presence of prominent emission features such as the Polycyclic Aromatic 
Hydrocarbons is clearly required to explain the observed colors for more than half of the sample. As a result, a general good 
agreement is observed between the data and predictions from the local starburst-dominated SEDs that have been used so far 
to constrain IR galaxy evolution. This is consistent with the star-forming nature of LIRGs derived from previous works, even 
though our approach cannot rule out the dominance of an AGN in some cases. Our study also supports the possibility of tracing 
the total IR luminosity of distant galaxies (up to z ~ 1) from their MIR emission. 
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1. Introduction 

The excess of faint sources in the deep extragalactic surveys 
performed in the MIR to sub-millimeter with their associated 
redshift distributions when available, together with the shape 
and intensity of the cosmic IR background (CIRB, Puget et al. 
1996), which measures the extragalactic IR light radiated above 
3 fim and is at least equal to the UV to near-IR extragalactic 
background light (Gispert et al. 2000; Hauser & Dwek 2001), 
suggest that a large fraction of the UV radiation of young stars 
was reprocessed by dust in the mid to far IR range over the 
Hubble time. The cosmic star formation history- the star for- 
mation rate (SFR) per unit comoving volume as a function of 
redshift- would be dominated by intense star formation phases, 
during which the bulk of the UV light is reprocessed by dust in 
the IR with a rapid decline of the SFR density since z ~ 1 and 
a flat or possibly slow increase from z ~ 5 to 1 (Chary & Elbaz 
2001, Lagache et al. 2004 and references therein). These mod- 
els called "backward evolution models" because they evolve 
the local luminosity function in the MIR as a function of red- 
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shift, assume that the shapes of distant SEDs in the MIR remain 
similar to the ones observed in local galaxies and are directly 
related to the total IR luminosity. 

The goal of the present study is to check the validity of 
this assumption up to z ~ 1 by combining data from ISOCAM 
(Cesarsky et al. 1996a), the MIR camera onboard the Infrared 
Space Observatory (ISO, Kessler et al. 1996), and MIPS, the 
Multiband Imaging Photometer for Spitzer (Rieke et al. 2004) 
onboard the Spitzer Space Telescope (Werner et al. 2004). 
Studying the MIR shape of distant galaxies is particularly rel- 
evant since direct FIR observations are rapidly confusion lim- 
ited and strongly affected by moderate sensitivity limits. On the 
other side of the peak of the IR SED, sub-millimeter observa- 
tions are presently limited to the detection of ultra-luminous 
IR galaxies (ULIRGs) above z ~ 2 (Chapman et al. 2003). 
In the near future, before the launch of Herschel and the ad- 
vent of ALMA, the best constraints on the evolution of lumi- 
nous IR galaxies (LIRGs), i.e. galaxies with Lir > 10 11 L Q , 
as a function of redshift and on their role in the global evolu- 
tion of galaxies will exclusively come from the MIR. Around 
z = 0.7, LIRGs are too faint to be measured through direct 
spectroscopy in the MIR with the infrared spectrograph (IRS, 
Houck et al. 2004a) onboard Spitzer. The only remaining tech- 
nique to contrain the MIR shape of distant LIRGs is therefore 
to combine MIR measurements for LIRGs located at different 
redshifts, hence spanning a larger range of rest-frame wave- 
lengths. In the present study, a redshift range of z =0.1-1.2 
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allowed us to study the 5 to 25 pm rest-frame part of the MIR 
SEDs. In this wavelength range (see Genzel & Cesarsky 2000), 
the SED is dominated by the combination of broad emission 
lines, generally interpreted as due to polycyclic aromatic hy- 
drocarbons (PAHs), and of the continuum emission of stochas- 
tically heated "very small grains" of dust transiently heated to 
temperatures of the order of 200 K. The thermal emission of 
"big grains" of dust heated to ~ 40 K also contributes partly 
to the MIR emission but peaks in the FIR (between 80 and 
100 pm typically) and contains the bulk of the luminosity radi- 
ated by galaxies above 3 //m. Finally, hot dust emission due to 
dust heated by an active galaxy nucleus (AGN) can also con- 
tribute to and sometimes even dominate the MIR spectrum of 
a galaxy. Before deriving the IR luminosity for a galaxy that 
will be used to compute its SFR, one must start by making 
sure that its MIR SED is not polluted by the AGN emission. 
The presence of PAHs as well as a rapid decline of the contin- 
uum emission below 5 pm strongly suggests a star formation 
origin for the emission (Genzel & Cesarsky 2000, and refer- 
ences therein). Local LIRGs are rarely affected by an AGN, 
contrary to ULIRGs above 10 12 3 L Q . Using a combination of 
template SEDs and deep X-ray surveys with XMM-Newton 
and Chandra, Fadda et al. (2002) derived an upper limit to the 
contribution of an AGN to the MIR light radiated by ISOCAM 
selected LIRGs up to z ~ 1 of 20 %. 

The MIR emission of local star forming galaxies- at 6.75, 
12 and 15 pm- was proven to correlate, with some scatter, with 
the total IR one, i.e. from 8 to 1000 /mi, which is largely domi- 
nated by the FIR component (Chary & Elbaz 2001, Elbaz et al. 
2002). However, the validity of these correlations in the more 
distant universe has not yet been established. Galaxies were 
less metal rich in the past. Their shape and compacity evolving 
with time might also affect their SEDs. Local galaxies SEDs do 
present some variations in the MIR as a function of metallic - 
ity. The metal-deficient (Z=Z /41) blue compact dwarf galaxy 
SBS 0335-052, for example, is bright in the MIR range but 
still does not show any sign of the presence of the PAH fea- 
tures which the authors interpret as an effect of the destruction 
of their carriers by the very high UV energy density (Thuan, 
Sauvage & Madden 1999, Houck et al. 2004b). 

Throughout this paper, we will assume H a = 70 km s _1 
Mpc" 1 , Q mat ter= 0.3 and Q A = 0.7. 

2. Sample selection and data reduction 

A 30'x30' field was imaged by ISOCAM at 15jt/m in 
the Marano FIRBACK region centered on (3hl3ml0s,- 
55°03'39"), about two thirds of which is covered by a larger 
MIPS-24;um image of 1.4° x 0.7° with a 80% complete- 
ness limit of 170//Jy. The central part of this field (7'x7\ 
UDSF for "ultra-deep survey firback") was covered down 
to an 80% completeness limit of 140//Jy at 15/zm with 
ISOCAM, a depth similar to the MIPS one (see Fig. [TJ. The 
ISOCAM- 1 5 jt/m UDSF belongs to the ISOCAM Guaranteed 
Time Extragalactic Surveys (IGTES, Elbaz et al. 1999). Six 
overlapping mosaics were performed in the micro-scanning 
mode of ISOCAM for a total of 2.7 hours of integration per 
sky pixel and leading to a final projected map with 2" pix- 




Fig. 1. Left: ISOCAM 15 pm image of the Marano FIRBACK 
field (central part of the ultra-deep survey, 7'x7') down to an 
80 % completeness limit of 140 piy. Black X represent the po- 
sitions of the common sources with MIPS plotted in the Fig. 13 
Right: MIPS 24 pm image of the same field down to an 80 % 
completeness limit of 170//Jy. White X represent the position 
of ISOCAM sources. 



els. The data were reduced using PRETI (Starck et al. 1996) 
and the transient correction as well as the assessment of the 
completeness limit were computed by Monte Carlo simulations 
with fake sources of known flux densities in the real images. 

The 24 pm observations were performed during the MIPS 
commissioning phase (IOC/SV) in November 2003. The scan 
map AOT was used, with an half-array overlap to cover about 
1300 Sq. Arcmin with high redundancy (20) and to get an 
integration time per sky pixel of about 230s (Papovich et 
al. 2004). The data were reduced using the Spitzer Science 
Center Pipeline and the BCD products (Basic Calibrated 
Data, Pipeline version S 10.0.3) were downloaded from the 
Spitzer archive 1 . PSF-fitting photometry was performed using 
DAOPHOT (Stetson 1987) with IRAF 2 . 

For the present study, we used a sample of 16 galaxies de- 
tected with both ISOCAM and MIPS at 15 and 24 pm respec- 
tively and identified with crosses in Fig. ^ This set of galax- 
ies was extracted from a larger sample of ISOCAM selected 
galaxies for which medium resolution (R=1200) VLT-FORS2 
spectra were obtained (Liang et al. 2004) in the "Marano 
FIRBACK" field. Oxygen abundances were obtained for half 
of them with values of 12+log(0/H) ranging from 8.4 to 9. 
With a median stellar mass of 5xl0 10 M and redshift of z ~ 
0.7, they are among the most massive galaxies in this redshift 
range but exhibit metallicities about twice smaller than their 
nearby counterparts suggesting that they will produce about 
half of their metals from z =1 to (Liang et al. 2004). Finally, 
it must be noted that even though no direct FIR measurement 
exists for this sample of galaxies at present, there is indirect ev- 
idence that the MIR-FIR correlations remain valid up to z ~ 1 . 
Indeed, MIR and radio measurements provide consistent pre- 
dictions for the FIR luminosity of distant galaxies up to z ~ 1 
(Elbaz et al. 2002, Appleton et al. 2004). 

1 http://archive.spitzer.caltech.edu, PID: 718 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation 
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Fig. 2. All SEDs sorted (from top to bottom) by decreasing red- 
shift. The SEDs are shifted by an arbitrary offset in vL v for vis- 
ibility; wavelengths are rest-frame. ISOCAM 15 /mi and MIPS 
24 //m luminosities are reported with the filter bandwidth and 
the 1-cr uncertainty. The label indicates the logarithm of the IR 
luminosity as well as the redshift and ID of each source, e.g. 
ID #4 is the galaxy UDSF04 in Liang et al. (2004). The filled 
circles are the luminosities that would be measured in both fil- 
ters for the plotted SED from the library of template SEDs of 
Chary & Elbaz (2001). Bold dashed line on gal.#31: SED of 
the AGN NGC 1068 normalized to best fit the measured 15 
and 24 /mi luminosities of the galaxy. Dashed line on gal.#l: 
SED of NGC 77 14 (Brandl et al. 2004). 

3. Discussion 

Fig-Elpresents the 15 and 24 fim luminosities of the 16 galax- 
ies detected with ISOCAM and MIPS in the UDSF fitted by 
a library of template SEDs from Chary & Elbaz (2001). This 
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Fig. 3. The observed colors vL v (15)/vL v (24) derived from 
our sample (open circles and vertical error bars, with num- 
bers refering to IDs from Fig. |2j, compared with the pre- 
dictions for different object types such as the Seyfert 2 
prototype NGC 1068, two enshrouded IR-luminous systems 
with deep silicate absorption (Arp220 and NGC 44 18), M82 
(log(LiR/L o )=10.6) and three template SEDs from Chary 
& Elbaz (2001, plain lines; the lowest line corresponds to 
log(L IR /L )=12.4). 

library of 100 template SEDs, from log(L IR /L )= 8.5 to 13.5, 
was generated between 0. 1 and 1000 fim to reproduce the trend 
observed for local galaxies between MIR and FIR luminosi- 
ties. The MIR (4-20 jum) part was produced by interpolating 
between the ISOCAM CVF spectra of four prototypical galax- 
ies (Arp220, NGC6090, M82 and M51), hence all SEDs con- 
tain PAH features. Because the 6.75, 12 and 15/im luminosi- 
ties of a local galaxy correlate with its total IR luminosity, Lir, 
it is possible to use a single MIR measurement to derive Lir, 
and therefore a SFR, for nearby galaxies. In "backward evo- 
lution models", such as Chary & Elbaz (2001), the 15 fim lu- 
minosity function is evolved as a function of redshift to repro- 
duce the observe galaxy counts and the CIRB, assuming that 
a template SED is uniquely determined by a given 15 fim lu- 
minosity. To fill the gaps between the observed 6.75, 12 and 
15 /mi correlations, the templates interpolated for all total IR 
luminosities were used. We used the same technique to derive 
Lj 1 ^ from L15 (measured with ISOCAM), Ljj from L24 (mea- 
sured with MIPS) and L^ 24 from the combination of L15 and 
L24. First, the 100 template SEDs were redshifted at the redshift 
of a given galaxy among the 16 and a flux density was deter- 
mined at 15 fim using the filter response of ISOCAM. Second, 
the SED which 15 pm flux density was the closest to the ob- 
served one was selected and normalized to reach the observed 
15 /mi flux density. The total IR luminosity of this galaxy is 
Lj 1 ^. Using the same strategy, we computed Ljg. Lastly, we se- 
lected the SED which minimized the^ 2 of the observed 15 and 
24 fim flux densities with their associated error bars and derived 
a third luminosity, Lj^~ 24 . The MIR luminosities in the rest- 
frame of the galaxies (crosses proportional to the filter width 
and error bar on the flux density) are compared to the SEDs 
used to compute L}^~ 24 in Fig.|5]and their associated MIR lu- 
minosities (filled circles). 
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The values of L IR derived by the three techniques are equal 
with a 1-cr dispersion of only 20%, which confirms the ro- 
bustness of this technique. Hence on average, the combination 
of both instruments will affect the predicted Lir values with a 
20 % dispersion only with respect to the use of only one of the 
two instruments. Hence the global shape of the template SEDs 
is such that it can be used to predict with this accuracy the other 
MIR value, suggesting that the SEDs did not vary very much 
since z ~ 1. In the "starburst" regime, i.e. below Lir = 10 11 L Q , 
the template SEDs do not provide a good fit to the galaxy 
(z = 0.4656, log[LiR/L ] = lO.7), which behaves very similarly 
to NGC 7714 (Brandl et al. 2004) as shown in dashed line in 
Fig. 13 The dispersion observed even among local galaxies, for 
example in the MIR-FIR correlations already suggested that an 
ideal library of template SEDs would have to include a varia- 
tion of shapes for each Lir bin. 

The bold SED in the middle of Fig.[2]presents a strong evi- 
dence for the presence of the 7.7 /urn PAH feature, band or com- 
plex. The rest-frame 9 fim luminosity of this 10 11 1 L galaxy 
(z ~ 0.7) is 1.8 times larger than the 14/im one, which is natu- 
ral when PAHs are present but rules out a hot dust continuum 
emission as the one locally observed in individual HII regions 
(such as M17, Cesarsky et al. 1996b), dwarf galaxies (Thuan, 
Sauvage & Madden 1999) or even Compton-thick Seyfert 2's 
(such as NGC 1068, bold dashed line in the Fig.fJ Le Floc'h 
et al. 2001). The SEDs were sorted as a function of increas- 
ing redshift in the Fig. |2] in order to illustrate the effect of 
the k-correction which acts as a low resolution spectrograph 
by shifting the broadband ISOCAM and MIPS filters to lower 
wavelengths with increasing redshift. The good quality of the 
fit was obtained without allowing the luminosity of the tem- 
plate SEDs to vary which illustrates the very good agreement 
with observations at 5-25 /urn rest-frame. 

The templates from Chary & Elbaz (2001) are therefore 
consistent with the observed 24/15 fim colors of distant sources 
considered as a whole. However it is worth mentioning that 
they do not provide a unique solution for every galaxy of the 
sample, since some of them can also be fitted by a feature-less 
continuum (e.g., ID #25 in Fig. 2). This issue is explored with 
more details in Fig. [3] (similar to Fig.l of Charmandaris et al. 
2004) where we compare the colors derived from our sources 
with what would be observed as a function of redshift for differ- 
ent object types such as the Seyfert 2 prototype NGC 1068, two 
enshrouded IR-luminous systems with deep silicate absorption 
(Arp220: Elbaz et al. 2002; NGC 4418: Spoon et al. 2001), 
M82 and three template SEDs (as in Fig|3). The comparison 
shows that galaxies with vL v (15)/vL v (24) < 0.8 can indeed be 
explained by power-law spectra with no significant contribu- 
tion from PAHs. This fitting degeneracy results in an additional 
uncertainty in the extrapolation to the total IR luminosity (typ- 
ically a factor of 2), and it also indicates that our approach is 
not sufficient in itself to discreminate between starbursts and 
AGNs. Nonetheless, a non negligible fraction of sources at 
z~ 0.2-0.8 (-40-50%) exhibit vL v (15)/vL v (24) colors larger 
than ~ 1 . The latter can not be reproduced by SEDs with- 
out strong emission features redshited in the ISOCAM band 
and boosting the flux at 15 fim. As predicted by the starburst- 



dominated templates, these colors are therefore the telltale sig- 
nature for the presence of PAHs in the distant Universe. 

4. Conclusion 

We have shown that by combining ISO and Spitzer in the pho- 
tometric mode, it was possible to constrain the MIR SED of 
distant galaxies too faint to be subject to direct spectroscopy 
in the MIR. Deep images with both instruments and at similar 
depths (140 and 170yuJy at 15 and 24/vm) do detect the same 
objects which illustrates the robustness of both instruments. 
The combination of 15 and 24 /mi flux densities measured with 
ISOCAM and MIPS was used to constrain the rest-frame 5- 
25 fim part of the SED of a sample of 16 LIRGs with z ~ 0.1- 
1 .2, taking advantage of the k-correction. Similar studies with 
larger statistics will become feasible soon by combining MIPS 
24/im with IRS 16yum peak-up imaging. Even though some 
sources were found to be either consistent with a starburst or 
an AGN-dominated SED, a significant fraction of our sample 
shows clear evidence for the presence of the broad bump asso- 
ciated to the 7.7 fim PAH feature in emission and for the silicate 
feature in absorption centered at 9.7 fim. Local template SEDs 
fitting the correlations between MIR and FIR luminosities pro- 
vide an overall good fit to the distant LIRGs supporting the pos- 
sibility to use them in models fitting galaxy counts and deriving 
a cosmic star formation history using a "backward evolution" 
of the MIR luminosity function. 
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